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Oh' The current voltage (IV) characteristics for perpendicular transport through two se- 

quentially coupled wires of finite length is calculated analytically. The transport within a 
Coulomb blockade step is assisted by plasmon resonances that appear as steps in the IV 
' characteristics with positions and heights depending on inter- and intrawire interactions. 

CNj , In particular, due to the interwire interactions, the peak positions shift to lower voltages 

' in comparison to the noninteracting wires which reflects the facilitation of current by in- 

' teractions. The interwire interactions are also found to enhance the thermally activated 

' ^ , current. 

fH ' Recent technological progress in creation of clean one-dimensional wires, either as carbon nanotubes^ 
, or as a part of semiconducting structures^ has attracted much attention to the study of transport through 
' quasi-one-dimensional systems, both experimentally and theoretically.^ ^ Transport through clean quasi- 
I one dimensional systems probes the main consequences of the Luttinger liquid (LL) theory®, the spin- 

^ ■ charge separation and the charge fractionalization. Experimentally, quasi- Id systems have a finite length, 
I or they are effectively divided into finite size fragments by impurities^. If the length is not very large, 

a the quasi-ld systems actually represent one-dimensional quantum dots, with finite charging energy and 
I ' discrete spectrum. Therefore, the theoretical study of finite-size effects on transport through wires of 
'"^ , finite length L that form the LL in the limit L ^ oo is of much importance for the interpretation of 
^ ■ experiments'". 

O ' In this Letter, we address those features of a LL that can be seen in experiments on perpendicular 
transport through two coupled quantum wires of finite length L, analyzing the current-voltage (IV) 
characteristics of the circuit consisting of the left reservoir, two quantum wires, and the right reservoir, 
coupled in sequence (see the inset in Fig. 1). The reservoirs are assumed to be noninteracting. There is a 
magnetic field applied perpendicular to the plane of the wires. Such system can be realized in form of two 
carbon nanotubes. A similar geometry has been employed in experiments on semiconductor structures^ 
with the difference that only one of the wires is of finite length there. Since the discreteness of the 
' spectrum is explicitely taken into account, the results presented here may find analogies in experiments 
on double quantum dots^". We do not consider spin effects, the Luttinger liquids are assumed to be 
P\| spinless. We derive an analytical expression for the current-voltage characteristics in lowest order in 
. tunneling to the reservoirs, concentrating on the region within a single Coulomb blockade step. 

We found that the interwire interactions facilitate transport at finite bias voltages, as well as thermally 
activated transport. Our results on IV characteristics are summarized in Fig. 1. The discreteness of the 
y spectrum of a finite-length wire transforms the power-law singularities, typical for the LL, into a series of 
^ [ discrete steps with envelope described by the power-law function of the infinite system^^. Each step in the 
(— ( IV characteristics corresponds to a plasmon resonance between the LL wires. The positions of the steps 
Q of the IV curve reflect the quantization of the charge-mode of the LL, in particular, the energies of two 
O . independent bosonic modes, propagating on the length L with velocities ui and U2.^ One of the velocities 
^ ■ (w2 in our case) always renormalizes down by interwire interactions, which results in the lowering of 
• i-H , plasmon excitation energies, and hence the facilitation of transport through excited states as compared 
' to the noninteracting wires. The height of each step reflects the strength of the corresponding plasmon 
^ ] resonance, and is also controlled by LL interaction constants. Therefore, IV characteristics should clearly 
■ - - ■ demonstrate the effect of LL interactions. The outlined features are rather insensitive to the geometry 
of the reservoirs. They can be seen in the case of two- and three-dimensional reservoirs contacted along 
a whole wire, as well as for a circuit with point-like coupling between the reservoirs and the wires, such 
as that realized in STM experiments^''^. In contrast, due to strong restrictions imposed by momentum 
conservation by tunneling, the interaction effects can hardly be seen in the case of one-dimensional 
reservoirs parallel to the wires. 

Our starting point is the Hamiltonian of two interacting Luttinger liquids of flnite length, with different 
Fermi-velocities vi and V2 and with a single-particle (IP) tunneling of amplitude t between the LL's. The 
length of the wires L is assumed to be large, so that the wave number k is conserved in the tunneling 
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process. In the presence of a magnetic field B perpendicular to the plane, the states of the wire 1 and the 
left reservoir acquire an additional phase factor e*^^ (b = eB/h, x denotes a coordinate along the wire) 
with respect to the states of the wire 2 and the right reservoir®. The application of the LL approach to 
wires of finite length is limited by the condition that the screening length for the Coulomb interaction 
in the wire is much less than the length of the wire L. This condition is fulfilled, for example, in recent 
experiments on semiconductors^'^. 

We take into account the electrostatic charging energy phenomenologically, shifting the groundstate 
energy by the capacitive energy E{N\,N2), where Ni,N2 denote the total number of electrons in the 
first and the second wire respectively. The calculation of E(Ni, N2) has been performed previously. ^■"'^"'^ 
Furthermore, it is assumed that the application of voltage between the wires shifts the energy of the 
one-particle states of the two wires by the amount eV relative to each other (band-shifting case^). 

Since transport between the LL's is determined by plasmon resonances, and in resonance the 
interwire tunneling has to be treated exactly, we first diagonalize the IP spectrum. The spec- 
trum of the diagonalized IP Hamiltonian consists of two branches with dispersions e±{k) = 

eV + {vi +V2)k± a/ {eV + {vi — V2)kY + /2. In general, the resulting spectrum is not linear any 

more for Vi 7^ V2. However, as long as |ey(ui — V2)k\ <C eV^ + Af^ one can expand the square root, and 
retaining only the first term of expansion consider the spectrum as approximately linear. Further we 
work in that approximation. Upon inclusion of two-particle interactions, the two branches e±{k) form 
two new interacting Luttinger liquids. 

After the diagonalization of the IP part of the two- wire Hamiltonian, the physical picture changes from 
the sequential circuit, the left reservoir - two coupled wires the right reservoir, to the circuit, where 
there are two LL's coupled in parallel to the left and right reservoirs with no direct tunneling between 
the Luttinger liquids. 

Note that two particle hopping terms between the new LL's in the particle-particle (pp) and particle- 
hole {ph) channels^^ are generated by the diagonalization of IP Hamiltonian. These terms will be omitted 
in further consideration. Their amplitudes are small at small tunneling t, though they play a crucial role 
close to the instability in the ph or pp channel (see below). 

After the IP diagonalization, the Hamiltonian of interacting LL's without the one-particle tunneling is 
of the form H = X^^^^ 2 Hi,+Hi2, where Hi, is a standard LL Hamiltonian without backward scattering^®, 

and the interwire interaction is of the form H12 = vi /_^^/2 If ■ ["UP1LP2R + HP1LP2L -\- L ^ R] : . ^2 
and 74 are the interwire interaction constants, j^{x) denote the chiral electron density operators of the 
two wires. 

In the bosonized form, the Hamiltonian looks like (here $ = ($1, $2) , © = (©i, ©2) ) 



L j1 



vigi ui(74-72)\ £r _ ( vi/gi ^1(74 + 72) 



where Ms — \ , x , Mg = \ . ' x , M , is the charge mode 

V^lW4-72) W252 / VWl(74 + 72) V2/g2 J 

velocity in wire v and is the LL interaction constant for wire v. Diagonalization of the Hamiltonian 
(1) is achieved by applying a composition of unitary rotations and subsequent rescalings, whereby the 
rescalings preserve the duality of the fields 4>i and 8^: $i ^ ^i^/TQ, Qi Qi/y/jQ. The same has been 
achieved by other means in®'^^'^^. The diagonalization renders (1) to the Hamiltonian of two uncoupled 
LL's with different charge mode velocities, ui and U2- H = X^^^^ 2^ / if {• i^xTliY '■ + • (^x^i)^ •} > 
where the dual fields rji, are related to the original ones by linear transformations $ = Q(p{VijV2)'^ , 
= Q6i(Ci, ^2)"^- In the case of identical wires, vi = V2 = v, 51 = 52 = g, the matrices (5(/>,e are given by 



i/g+74-1-72 



9+74 — f2 
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1/g — H — f2 



Ql' = 72 ^^f^ ' Qf - -Ql' - 775 -ir^ ^ and = 1/(2Q;^). The behavior of the 
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ff-74+72 

renormalized velocities Ui and U2 versus the interwire coupling strength at different values of intrawire 
interaction constant is shown in the inset in Fig. 1. Each velocity determines the inter level spacing 
for a branch of the plasmon spectrum by the relation Aj = TrUj/L, (i = 1,2). The velocity U2 that 
renormalizes down by interwire interactions corresponds to the mode, where the charges of two wires 
fluctuate in antiphase (the pseudospin mode in the case of identical wircs®'^^'^®). There is a special point 
in the dependence of the renormalized velocities on interwire interactions, where ui = U2- For identical 
wires, this point is given by the relation (74 + 72)51 = (72 — 74)/5- The farther away from that point, 
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the lower is the velocity U2- Moreover, U2 turns to zero when one of the matrices, or Me becomes 
degenerate. The latter happens when interwire and intrawire interactions are comparable. Particularly, 
in the case of identical wires, and for 72 = 0, U2 = is reached for 74 = vg. The interwire tunneling 
in coupled infinite LL's is known to result in an instability either in the ph or in the pp channel at 
strong enough interwire interaction^^. We expect that one of those two instabilities occurs at the point 
U2 = 0, even though the length of the wires is finite. Indeed, the velocity U2 changes the effective length 
of the wire for the corresponding plasmon mode by the obvious relation L Lv/u2, hence the physics 
of an infinite LL should be recovered for U2 — > 0. Complete analysis of the instability employing the 
renormalization group treatment is left for future publication. 

Now we turn to the implications of the above results for the current- voltage characteristics. We consider 
here the case, when tunneling between a wire and the next reservoir is much larger then tunneling between 
the wires. Then the whole voltage drop takes place between the wires. Increasing the bias voltage, the 
energy difference between the two IP spectral branches e+{k) — t-{k) increases, which results in the 
suppression of tunneling matrix elements Ti^r and T2^i (i = 1, 2 relates to the LL channel, v = r,l relates 
to the reservoir), and hence to the suppression of current. However, for a given chirality X (x = 1 for the 
left, X = ~1 for the right chirality), that effect can be compensated by changing the external magnetic 
field according to b = —xeV/vi (then both LL channels have the electrochemical potential c;losc! to the 
left reservoir) or 6 = xe V/ V2 (then both LL channels have the electrochemical potential close to the right 
reservoir). In what follows we consider the case b = —eV/vi. 

The single-particle tunneling current / under the applied voltage V is given by the sum of the currents 
through the two parallel LL channels. In each channel, there are contributions from states with the 
left and the right chiralities. The calculation of the current is outlined in Ref. [^^], whereas we use 
the bosonization formulas for a finite length wire with free boundary conditions to calculate the fermion 
correlation function in each channel. Furthermore, we describe the reservoirs in terms of their density 
of states with energy e, Pr,i {c)- That is appropriate in the case of two- and three-dimensional reservoirs 
with large bandwidth, where the momentum conservation by tunneling does not impose a restriction on 
energy. The current-voltage characteristics is obtained as follows. 



where 6,1 = (Qf)^ + {Qf)\ 6,2 = {Q'^f + (Of )^ = E{N, ± 1, TV,) - E{N,, N,) ± n(v, + nu, + 



mu2)/L — eV, and a is a short-range cutoff. Each term in the sum over n, m describes the contribution 
from a plasmon resonance between the two wires to the current. The tunneling density of states represents 
a superposition of two sequences of equidistant (^-peaks with periods nui/ L and 1:112/ L. Both the position 
and the weight of each peak is completely defined by interactions. The structure of the tunneling density 
of states has much in common with the structure of photoluminescence spectra from highly excited 
nanorings, described in detail in Ref. \^^]. Indeed, in both cases the correlations of two finite size LL's 
determine the relevant physics with the difference, that it is the four-point correlation function that 
determines the luminescence spectra. 

A set of IV characteristics for different values of interwire interaction 74 is shown in Fig. 1. The 
shift of the position of the second step to lower voltages with increasing 74, which corresponds to the 
current through the first excited plasmon state, demonstrates the facilitation of current by interwire 
interactions. For single wall carbon nanotubes, the charging energy is 3 to 6 times the noninteracting 
interlevel spacing, ^'^'^ independently of the nanotube length. Therefore, the plasmon resonances with 
m, n ^ 3 — 6 arc expected to be seen experimentally within a single Coulomb blockade step 

For one-dimensional reservoirs, the momentum conservation by tunneling fixes the energy difference 
between the modes in the reservoir and in the two-wire system. The IV curve yields a sequence of peaks. 
The positions of the peaks are determined by interference phenomena, which are only partially controlled 
by interactions. Moreover, even applying an external magnetic field, it is generally impossible to satisfy 
the momentum conscirvation and to avoid the diminishing of tunneling by the bias voltage simultaneously. 
We conclude that the interaction effects can hardly be seen if the reservoirs are one-dimensional. 

Another effect of interactions is in facilitating plasmon assisted transport at small voltages and finite 
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temperatures. For small temperatures, we obtain 7(T) = Iir{T) + hiiT), where 



l^|7i,.|V(eF)exp(-ff; 



(3) 



The enhancement of the thermally activated current can be seen by comparison of Eq. (3) with the result 
for the system without the interwire interactions that reads 7,^ oc exp (^7?) / [l — exp (^)] • Due to 
the smaller interlevel distance for plasmon excitations in the interacting system (smaller velocity ^2), the 
thermally activated current is enhanced. 

As the singularity U2 ^ is approached, the distance between the m and m+1 step in zero-temperature 
IV characteristics diminishes, the height of each step grows and formally diverges with 6^2 at the point 
U2 = 0. The thermally activated current also grows and diverges as a power-law (T/w2)'''^ . To describe the 
behavior close to the singularity, the ph and pp hopping operators need to be included in the theoretical 
treatment, which is beyond the scope of the present paper. Beyond the transition, the fluctuations of the 
fields (772, ^2) are frozen and do not contribute to the current. The charge fluctuations actually correspond 
to a single LL of finite length. The IV characteristics is given by Eq. (2) with U2 = 0, 6^2 = 0. The 
thermally activated current is suppressed compared to the noninteracting system, since ui > vp. 

The author acknowledges illuminating discussions with M.E. Raikh, who initiated this work. Discus- 
sions with D. Pfannkuche, H. Heyszenau, and S. Kettemann are gratefully acknowledged. 
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FIG. 1. IV characteristics for difFcrcrit intcrwirc interaction constants 74. The position of the second step 
reflecting the current through the first excited plasmon mode shifts to lower voltages with growing 74, which 
indicates the facilitation of current by interactions. In the shown voltage range, the curves with 74 = 0.2; 0.3 show 
equidistant steps, reflecting the contribution of the mode U2 only, whereas both modes ui and U2 contribute to 
the curve with 74 = 0.1. For 74 = 0.5, iti ~ U2, the positions of the stops from both modes coincide, the height 
of each step is the sum of the two contributions. Normalization of the current Iq = lQixe\Ti^r\^ pr^Ep) / a. Inset: 
the velocities of plasmon modes wi and W2 vs 74. Other parameters: g\ = 0.28, 52 = 0.26, 72 = 0.05. 
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